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Exact and asymptotic lineshape expressions are derived from the semi-classical density matrix 
representation describing a set of closed three-level A atomic or molecular states including deco- 
herences, relaxation rates and light-shifts. An accurate analysis of the exact steady-state Dark 
Resonance profile describing the Autler-Townes doublet, the Electromagnetically Induced Trans- 
parency or Coherent Population Trapping resonance and the Fano-Feshbach lineshape, leads to the 
linewidth expression of the two-photon Raman transition and frequency-shifts associated to the 
clock transition. From an adiabatic analysis of the dynamical Optical Bloch Equations in the weak 
field limit, a pumping time required to efficiently trap a large number of atoms into a coherent su- 
perposition of long-lived states is established. For a highly asymmetrical configuration with different 
decay channels, a strong two-photon resonance based on a lower states population inversion is estab- 
lished when the driving continuous-wave laser fields are greatly unbalanced. When time separated 
resonant two-photon pulses are applied in the adiabatic pulsed regime for atomic or molecular clock 
engineering, where the first pulse is long enough to reach a coherent steady-state preparation and the 
second pulse is very short to avoid repumping into a new dark state. Dark Resonance fringes mixing 
continuous-wave lineshape properties and coherent Ramsey oscillations are created. Those fringes 
allow interrogation schemes bypassing the power broadening effect. Frequency-shifts affecting the 
central clock fringe computed from asymptotic profiles and related to Raman decoherence process, 
exhibit non- linear shapes with the three- level observable used for quantum measurement. We point 
out that different observables experience different shifts on the lower-state clock transition. 

PACS numbers: 32.70.Jz, 32.80.Qk, 37.10.Jk, 06.20.Jr 



I. INTRODUCTION 



In the 1930s, molecular-beam magnetic resonance 
techniques achieved very high precision allowing the 
observation of atomic/molecular systems essentially in 
total isolation [l[. The Rabi method revealed coupling 
interactions between internal energy states and provided 
plenty of information not only on atomic and molecular 
structure, but also on nuclear properties. In the 1950s, 
N.F. Ramsey realized a scheme with much higher 
resolution by increasing the interaction time between the 
atom or molecule and the oscillating field 0. Still today, 
this technique provides the highest resolution in order 
to follow a dynamical evolution of wave functions and 
probe its phase accumulations. Control and elimination 
of systematic frequency shifts dcphasing a wave function 
oscillation at a natural Bohr frequency are fundamental 
tasks to achieve precision measurement @ . 



An alternative tool to probe by the Rabi or Ram- 
sey sequences, for a dipole-forbidden transition, is to 
radiatively mix the atomic or molecular states. As an 
example, without natural state mixing from spin-orbit 
interaction, a long-lived Raman coherence between a 
ground state and a long-lived (as in alkalis) or metastable 
(as in two-electron atoms) level, to be referred as clock 
states, can be established by a two-photon process via an 
upper excited level, thus forming a three-level A system. 
The properties of such a system are strongly determined 
by an optical pumping mechanism leading to a formation 
of a Dark Resonance associated with a trapping of the 
atomic population in a coherent superposition of states 
fflioj. Since such quantum superposition states are 
radiatively stable, the associated Raman coherence 
production leads to extreixiely narrow Dark Resonances 
allowing high-resolution frequency measurements [11 
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Such coherences were explored for single trapped ions 
[l^ ITsj , microwave clocks []3 , microwave chips [l^ [l^ , 
optical lattice clocks [17 L multi-photon excitations 
[isl [l9| , or nuclear clocks l^j . Similar coherent super- 
positions are used in solid-state physics for quantum 
information [2l|, in super-conducting circuits [l^, l23| . 
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in a single impurity ion inserted into a crystal [241 ] . 
in quantum dots [25| with protection against random 
nuclear spin, interactions [26|, and in opto- mechanical 
systems [23, H^- They are also actively considered 
within the future challenge of realising nuclear systems 
for quantum optics in the X-ray region (29l . Isoj . Dark 
Resonances, largely exploited in Quantum Optics, 
have been extended to the pre par ation of molecules 
in ro- vibrational ground states [31[, and to coherent 
superposition of atomic-molecular states in order to 
measure atomic scattering lengths and lifetimes of exotic 
molecular states (32| - [35j . 

Three-level narrow resonances are associated to 
quantum interferences produced by amplitude scattering 
into different channels and are strongly dependent on 
the configuration intensities and detunings [sgI . [37| . 
Such different lineshapes were associated to quantum 
interferences as the Autler-Townes (AT) splitting of the 
resonance [s^, the Fano-Feshbach (FF) canonical form 
(39| , the Dark Resonance (DR) lineshapeH known as 
Coherent Population Trapping (CPT) [a, l4p| or Electro- 
magnctically Induced Transparency (KIT) [4ll. l42j . 

Accurate calculation of the lineshape for a quantum 
superposition resonance requires numerical integration of 
the Bloch's equations. The literature reports on efforts 
to establish approximate analytic equations applicable 
to each particular case and, in a few cases, exact but 
rather complex [l^, lisl - lisj . The application of A 
schemes to high accuracy atomic clocks, in microwave 
or optical domains, requires to determine precisely the 
physical processes affecting the resonance lineshape 
and the shifts of the clock frequency. For that purpose 
this work provides a careful analysis of the lineshape 
dependence on different parameters characterizing the 
atomic or molecular system under investigation. 

The standard clock interrogation of a three-level 
system involves continuous excitation of the two lower 
states while sweeping through the Raman resonance. For 
that regime, starting from the steady-state analytical 
solution of three-level Optical Bloch Equations wc derive 
the exact expression of the resonance lineshape where the 
role of the relaxations and dephasing rates determining 
the absorption profile is expressed with physical mean- 
ing. Our detailed discussion of key linesha pe p arameters 
expands previous analysis 0, [lfl,[ll,|3l,l33l- We show 
that Autler-Townes, CPT, KIT and FF lineshapes are 
associated to a universal two-photon resonance lineshape 
depending on system parameters [49j . The analytical 
expressions for the frequency shift associated either 
to the FF extrema or to the FIT resonance point out 
dependencies not obvious in a perturbation treatment. 



An alternative clock operation scheme is based on a 
Raman-Ramsey scheme with the application of time sep- 
arated but resonant two-photon pulses, experimentally 
introduced in the microwave domain [SOl |5l| a-nd ex- 
tended to rubidium cold atoms 



beam of sodium atoms [53| . While in the standard Ram- 
sey approach, a coherent superposition of clock states 
in the bare atomic/molecular basis is dynamically pro- 
duced by a 7r/2 pulse depending on pulse duration and 
laser power, the coherent superposition, in the three-level 
two-photon approach, is created by an optical pump- 
ing process long enough to reach a steady-state. This 
scheme overcomes the power broadening mechanism of 
the continuous wave resonance allowing to obtain high 
contrasted signals in a saturation regime. This idea was 
extended in refs. [s^ - iSGj to the realization of FIT- Raman 
(and hyper-Raman) optical clocks with alkaline-earth- 
metal atoms. The time-separated and individually tai- 
lored laser pulses may be designed to create an atomic 
coherent superposition while eliminating off-resonant ac 
Stark contributions from external levels modifying the 
optical clock resonance [53|. For the regime of the first 
laser pulse long enough to produce an efficient coherent 
superposition, we present here a detailed analysis describ- 
ing the dependence of the DR lineshape on the system 
parameters. 

Within the quantum clock framework, the determina- 
tion of lineshapes and resonance shifts in different ex- 
perimental configurations remains an important issue, to 
be carefully investigated within the present work. An 
important result of the present analytical and numerical 
analysis for the resonance frequency-shift of a three-level 
quantum clock, is that different lineshapes versus the op- 
tical detunings arc obtained depending on the experimen- 
tally detected population or coherence observable. 

The A system and the Bloch's equations for an ho- 
mogeneous medium are introduced in Sec. II, where an 
adiabatic analysis of the time dependent equations de- 
termines the approximated time scale required to pro- 
duce an efficient atomic/molecular coherent superposi- 
tion. Sec. Ill establishes an exact treatment of the ex- 
cited state steady-state regime and derives the key in- 
formations on the Dark Resonance lineshape. In Sect. 
IV, we derive the steady-state profile of a two-photon 
resonance between clock states. In Sec. V, we focus our 
attention to the Raman coherence lineshape observed be- 
tween clock states. Finally Sec. VI analyzes DR fringes 
produced with resonant two-photon pulses separated in 
time mixing steady-state properties and Ramsey oscilla- 
tions. A detailed analysis of the fringe properties is de- 
rived in the adiabatic regime where the first pulse estab- 
lishes a steady-state solution and the probe pulse dura- 
tion vanishes. Instead only dynamical properties of these 
phase-shifts were demonstrated in refs [HJ, [s^ [72|, [73| . 
In the Appendix A, we rewrite lineshape population so- 
lutions in terms of generalized multi-photon transition 
rates enhancing one and two photon transition rates in 
the three-level system. We finally derive in Appendix B a 
first order analytical expression of the central fringe Ra- 
man frequency-shift associated to the pulsed Dark Reso- 
nance lineshape. 
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inspired by the Ezckicl's group work at MIT on a thermal 



3 



7l /s/^i^i 

1^ 



FIG. 1: (Color online) Closed three-level A configuration in 
the density matrix representation including relaxation rates 
r3i,r32, and decoherences 71,72,7c- Optical detunings are 
Ai, A2. The parameter Sr — Ai — A2 defines the Raman 
resonance condition. Qi and 112 define the couplings with the 
applied laser fields. |1) and |2) are the clock states while |3) 
is the excited state. 



decay configuration. Optical coherences are relaxed with 
terms 71 , 72 where 71 -|- 72 = 7- In a pure radiative pro- 
cess [53] , optical decoherences are related to spontaneous 
emission rates by the relation 71 = 72 = {Tsi + r32)/2. 
The pi2 decoherence is described by the 7c dephasing 
term. The optical Bloch equations describe the temporal 
evolution of the density matrix elements in the RWA as 



dp 
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dt 

dp22 

dt 

dp33 

dt 
dpi3 
dt 

dp23 

dt 
dpi2 
I dt 



- 2f}i/m{pi3} + Tgipag, 

- 2n2lm{p23} + r32P33, 

217i/to{pi3} + 2172/m{p23} - ^P33, 

- [71 + iAi]pi3 + iCl2Pl2 - i^i{p33 - Pll), 

- [72 + iA2]p23 + 'i^lP21 - i^2iP33 - P22), 

- [7c + i(Ai - A2)]/9l2 + ifl2Pl3 - i^lP32, 



(4) 



II. THREE-LEVEL OPTICAL BLOCH 
EQUATIONS 

The Doppler-free three-level system presented in Fig.[T] 
is described by the density matrix pij — 1,2,3) obey- 
ing the Liouville equation 



dt 



Pij 



Pi,kHk,j 



) + np,. 



(1) 



The coupling of the three atomic or molecular states to 
two coherent radiation fields, see Fig. [1] is described 
within the rotating wave approximation (RWA) by the 
following Hamiltonian: 



^Ai ni 
H = h\ A2 ^2 

02 



(2) 



where Ai,A2 are the detunings of the two fields. De- 
pending on the transition the Rabi frequencies fli and 0,2 
driving the system are determined by the product either 
between electric dipole and electric field amplitude or be- 
tween magnetic dipole and magnetic field amplitude. It 
is worth noticing that Rabi frequencies defined here are 
half of the definition of ref. [131 • The matrice TZpij taking 
into account relaxation and decoherence phenomena, is 
written 



'r3iP33 -7cPi2 -7iPi3\ 

T^Pij = I -7c/321 r32P33 -72P23 

-71P31 -72P32 -r/933 J 



(3) 



The total spontaneous emission rate F is composed by the 
r3ii r32 rates (with F = F31 -I-F32) describing cither alka- 
line (F31 « F32) or akaline-earth (F31 F32) three- level 



with pji ~ p*j . The population conservation of the closed 
system is given by pu + P22 + P33 = 1- The optical 
detunings will be related to the Raman detuning dr by 
Ai = Aq and A2 = Aq — Sr- Aq represents the common 
optical detuning for a configuration where one laser is 
fixed while the other is frequency scanned. Notice that 
within the approach of deriving the Ai, A2 laser modes 
from a single source by modulation at frequency Sr, with 
Ai = Ao ± Sr/2 and A2 = Aq T Sr/2, the light-shift 
derivation should be modified. 

Eqs. ([4]) describe the transient dynamics and the 
steady state of populations and quantum coherences. 
A complete state mixing is reached when all atoms or 
molecules have been pumped efficiently into the dark 
state, a coherent superposition of the lower states. Thus, 
a pumping time is required to achieve an optimal atomic 
fraction trapped into the coherent state superposition. 
We derive such time scale evolution from Eq. (j4]) by an 
adiabatic elimination of the time derivative, dp33 (t) /dt = 
dpi3{t)/dt = dp23{t)/dt = for pulse durations greater 
than F;^;^^ , V^/^ i ^ the population ^33 and optical coher- 
ences P23 and pi3 evolve more quickly than the popula- 
tions pii, P22 and the Raman coherence pi2. 

We investigate the dynamics of the three-level A sys- 
tems using various combination of long and short two- 
photon pulses separated in time. A straightforward tem- 
poral analysis of the resulting adiabatic set, similar to 
NMR equations [s^, [s^, exhibits two damping times: 
Tosc, determining the phase memory of the Raman co- 
herence precession (equivalent to a transversal or spin- 
spin relaxation rate), and Tp which determines the typi- 
cal population transfer into the dark state superposition 
(similar to a longitudinal or spin-lattice relaxation rate) 
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[60|. At low optical saturation fii, ^ T, 7, we have 

-1 



Tosc{Aa) 
^p(Ao) 



) ~ 7c + 



71 



Of 

72 

h ^ 
72 



7i 

-1 r 



1 + T 



f^i/72 - f^?/7i 



_1 (5) 



where we have introduced the following generalized re- 
laxation rates: 



A?+7f 

7» 



1,2. 



The T generalized branching ratio difFcrcnce is 
f(r!?/7i-f^i/72) + T 



T 



m/ii 



1 



(6) 



(7) 



with the T normalized branching ratio given by T = 
(Fai -r32)/(r3i -l-r32). The Tosc and Tp timescales play 
a key role on the population transfer between atomic or 
molecular states [5J| . Indeed, optical coherences are effi- 
ciently generated only when the J7i and Rabi frequen- 
cies are applied for a time r exceeding , ■ ^ short 
pulse duration having V^-^ , < t < Tp will be Fourier 
limited and will lead to a weak contrast resonance profile, 
whereas a long pulse with r ^ Tp 3> T'^i , ^i^^ elim- 
inate all time dependencies in lincshape and frequency 
shifts. This regime is latter examined within the follow- 
ing section. 
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III. THE STEADY-STATE LINESHAPE OF pss 

A. The Dark Resonance 

In examining the steady-state situation with all time 
derivatives in Eq. (|4]) set to zero, we find the exact ex- 
pression for the population of the upper state P33 



P33 



(Ai 



A2)2 



■ 7c7e// 



(Ai-A2)[Ai-A2 



A/] + r^,/ 



(8) 



where Ai — A2 = (5r is the Raman detuning. The expres- 
sion of A/ , the frequency shift affecting the Raman de- 
tuning, and Fe//, halfwidth of the two-photon resonance, 
are reported in the following subsection. The coherence 
decay rate is 7e// = 7c + 7* with the 7* saturation rate 
of the Raman coherence given by 



1 = — + —, 

72 71 



(9) 



Notice that 7e// coincides with t^^^ for the Ai = A2 = 
resonant laser case. 

of Eq. ([5]) represents the signal amplitude. It con- 
tains the broad features of the P33 dependence on the 
optical dctunings Ai, A2. The fraction on Eq.®, whose 
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FIG. 2: (Color online) Three-level spectra vs 5,- Raman 
detuning observed on the pzz population. For all spectra 
Fsi = = r/2, 7c = 0. (a) AT spectrum at Ao = 0, 
Q.1 = 5r, and 0,2 = F. (b) Dark/EIT resonance at Ao = 0, 
fii = 0.2F, and ^2 = 5 x lO'^F. (c) FF resonance for Ao = F, 
Q,\ — 0.2F, 0,2 = 10~^r. Solid lines from the analytic solution 
Eq. ((8]) and dots from the numerical integration of Eq. Q at 
r > rp(0). 



values lie in the [0,1] interval, determines the P33 nar- 
row variation with the Raman detuning Ai — A2. The 
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three-level signal amplitude S is given by 

S 



8"^ = 



l + (3-E)5 + 



7172 



(10) 



where we introduced the normalized dimensionless Rabi 
frequencies [6l| 



r327i 



2 

^2 = ^27^ 



2 r3i72 



(11) 



f^?r3272 + ^jr^ai ' 



The saturation parameter 5* driving the population ex- 
change between energy levels is determined from Ein- 
stein's rate equations as: 



S 



(r3272f^? 



^3111^1)- 



(12) 



The imaginary parts of optical coherences are related to 
the excited state lineshape expression by the relation: 



Im{p^3} 



2a 



P33 (i = l,2). 



(13) 



Therefore their lineshape is equivalent to that of P33. 
In Appendix A, we recast all population lineshapcs in 
terms of multi-photon transitions rates, pointing out the 
light-shift contributions to the optical detuning terms. 
We verified that the numerical results to be presented in 
the following can be derived also from that solution. 

Depending on the detuning and intensity of the lasers 
or microwaves driving the three-level system, the line- 
shapes associated to Eqs. ([5]) and (|13p present very differ- 
ent features, known as AT spectra, Dark/EIT resonance, 
and FF profile, associated to different degrees of interfer- 
ence between two-photon transition amplitudes. 

The AT profile appears when one Rabi frequency is 
much larger than the natural linewidth of the excited 
state (fi; ^ r, i = 1,2), and in addition fti ^ rt2, or 
viccversa. Two splittcd resonances, a doublet structure, 
appears in the frequency spectrum, as shown in Fig. HJa) . 

At Rabi frequencies smaller than the excited state 
width, we reach the DR or EIT configuration where 
a narrow two-photon resonance is established from the 
quantum destructive interferences between the transition 
probability amplitudes [10,1111 as seen in Fig.[2][b). The 
system is placed in the dark state uncoupled from the 
driving fields. Note that at exact resonance, Ai = A2 = 
0; P33 — when 7c = 0. These regimes are characterized 
(Ai w A2 < r), with Oi « ri2 for the Dark Resonance, 
and rii ^ ilj; Hi, flj < T for the EIT resonance. 

The FF lineshape plotted in Fig. ^c) is originated 
when one Rabi frequency is much larger than the sec- 
ond one, and in presence of an optical detuning from the 
excited state(ili ^ ^Ij, Ai,A2 > F). Two resonances 
appear in the FF spectrum, one broad corresponding to 
the saturated one-photon resonance. The second sharp 



feature exhibits a characteristics asymmetric response, 
highly sensitive to changes in the system parameters and 
centered around the Sr ~ Raman detuning. Its mini- 
mum is associated to the DR, or EIT dip, while the max- 
imum is the Raman peak, or bright resonance, associated 
to the preparation of the coherent superposition of |1 > 
and 1 2 > states coupled to the driving electromagnetic 
fields. Both the EIT dip and the Raman peak are mani- 
festations of the interference between the one-photon and 
two-photon amplitudes [IBl. The asymmetry of the FF 
profile is reversed by changing the relative ratio between 
the Rabi frequencies. 



Fe/j Raman linewidth 



The sub-natural EIT resonance of Fig. ^h) experi- 
ences a linewidth which is power broadened by the op- 
tical saturation rate 7* of Eq. ([9]). Let's note that such 
power broadening is important even for a laser intensity 
where the saturation is negligible on the optical transi- 
tions. In fact ref. (62j introduced a coherence saturation 
intensity, defined by 7* = 7c, smaller than the optical 
saturation intensity. The exact expression of the Fe// 
Raman halfwidth is 



e// 



(7e// -I- 37c5') 7e// 



Xy/lTC, (14) 



where the factor C, is 



F^ [2AiA2 + A?-Al]+7; 



(71-12)8 /p 
27172 ^ 



ileff 
3272 



■ 37c5)7c// 
r3i7i) 



(15) 



7e// + 37cS' 



C is very small (C <C 1) in a quasi-resonant laser inter- 

7c 7c 

action (according to the condition — A^,— Aj <C i^^) 

and a pure radiative process. As long as 7c, Qi,Q2 F, 
Fe/ / determines the half- linewidth of the sub-natural res- 
onance. In that regime and for the pure radiative case, 
Fe// is well approximated by rJ^J,. 



C. A/ Raman frequency-shift 

The frequency shift A/ correcting the 6r — Raman 
detuning condition is given by 



A, 



^LS + 'AOS, 



(16) 
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FIG. 3: (Color online) Frequency shifts of the Dark/EIT dip and Raman peak observed on the P33 excited state population 
versus the optical detuning Ao from Eq. (|19[) (solid blue and red lines, respectively), and from the numerical integration of 
Eq. (Q (solid dots • and open dots o, respectively). Parameters are Fai — — r/2, fli — 5 x lO^^F, and Q2 = lO^'^F. In 
(a) ■jc — 5 X 10~*F; in (b) ■jc ~ leads to a vanishing shift for the EIT resonance. 



with the Als light-shift (LS) expression including the 
saturation effect given by 
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7172 



1'" ' 7172 



l + (3-^)5 
(71-72) r5'[72Ai 



- 71A2 



(17) 



271727 



l + (3-^)5- 



7172 



Notice that the second term of the above expression van- 
ishes for the symmetric A scheme, i.e., with 71 = 72 
(pure radiative process). In that case, the first term of 
Als could be associated to the light-shift expression, as 
pointed out in The Ads decoherence shift (DS) de- 
pends on the 7c rate as 



A 



DS = 7c 



s [r(Ai + A2)-(r3iAi -r32A2)] 



27172 



l + (3-^)S- 



7172 



(18) 

Let us emphasize that A/ is always null at Ai = A2 = 
resonance [6311. 



D. Approximated frequency-shifts of EIT and FF 
resonances 

Instead of the previous Section exact expression cor- 
recting the Raman detuning condition in the denom- 
inator, it is useful to derive the effective shift of the 
Dark/EIT resonance minimum, which is very relevant 
for precision spectroscopy or clock resonance. The calcu- 
lation of the DR/EIT and FF frequency shifts requires 
to examine Eq. ([5]) with the Raman detuning 8^. as a 
free parameter. A valid approximation for the EIT and 



Raman-peak shifts in various excitation configurations 
can be found when optical detuning Ai ~ A2 ^ Aq are 
tuned around the Raman condition Ai — A2 = 5r- A 
differentiation of Eq. (|5|) versus the 5^ parameter leads 
to roots of a quadratic equation defining the following 
extrema (Sw33(Ao) of the EIT/FF lineshapes: 



<^W33(Ao) ~ 




(19) 



The =F solutions refer to the extrema of the FF line- 
shape. The minus (plus) sign holds for EIT dip (Raman 
Peak) when jcjeff/^lff < 1, and the opposite when 
■^cjeff/^lff > 1- Fig. [3] shows the continuous-wave 
frequency shift versus the common optical detuning Aq 
using Eq. dH]). Results for both 7^ = S.IQ-^F and 7c = 
are presented. The EIT decoherence shift, proportional 
to 7c exhibits a linear dependence on the optical detuning 
when radiative decay terms are symmetrical. As pointed 
in p^ . when the decoherence term vanishes (7c = 0), 
there is no shift of the EIT minimum [g^l except if we 
take into account external small off-resonant level con- 
tributions [65| . The Raman peak shift produced by light 
and decoherence varies with the inverse of the optical de- 
tuning as discussed by [s^] and observed experimentally 



12| . For a quasi- resonant interaction we can further 



simplify the above expression for the Raman Peak and 
the EIT dip when 7c7c///r2^^ < 1 and A^/F^^^ < 1. 
Near the two-photon resonance, the shifts of Eq. (fTO)) can 
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be accurately approximated as 




-15 -10 -5 5 10 15 
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r2 

.RPf ^ \ r^O eff 



(20) 



The EIT dip frequency-shift is thus roughly given by the 
product of the 7c decoherence rate and the A/ Raman 
frequency shift divided by the Fg/ / linewidth of the sub- 
natural resonance. Such a dependence was pointed out 
by and by [1^ based on a theoretical analysis of , 
and as mentioned in l48i. was earlier derived in 16711 . 
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IV. STEADY-STATE LINESHAPES OF CLOCK 
STATE POPULATIONS pii,p22 

A. The two-photon resonance 

We focus now on clock-state resonances observed on 
the pii, P22 populations and linked to the Raman coher- 
ence between those states. The clock state populations 
may be expressed in an exact form similar to that of 
Eq. dH]) as 



Pii = S 



11 (Ai - A2)[Ai - A2 - All] 



P22 S 



22 (Ai - A2)[Ai 



A2 + A22] 



7?i 



(Ai-A2)[Ai-A2-A/] + r2^/ 



7I2 



(Ai-A2)[Ai-A2-A/] + r2^.^ 



(21) 



(c) 
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FIG. 4: (Color online) Three-level spectra versus the 5r Ra- 
man detuning observed on the pn population using Eq. (|21|l . 
solid blue line, and the numerical integration of Eq. solid 
dots (•). In (a) AT resonance, in (b) Lamb dip lineshape, 
and in (c) FF resonance. System parameters as in Fig. 2. In 
(a) and (b) Rabi frequencies as in Fig .2. In (c) f2i — 0.2T, 
n2 = 5x IQ-^r, and Ao = F. 
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FIG. 5: (Color online) Population inversion resonance be- 
tween the |1 > and |2 > states, monitored on pn and P22, for 
unbalanced decay rates, F32 = 10~^F, Fai = F — r32. Solid 
dots (•) are the result of the numerical integration of Eq. 
Other parameters are Ao = 0, 7c = 0, fii = 5 x 10~^r and 
fJa = 5 X 10-*F. 
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FIG. 6: (Color online) Frequency shift of the population inversion resonance observed on pii(or P22) as derived from Eq. (|25|) 
(solid blue line) and frequency-shift of the P33 EIT resonance from Eq. (|19p (solid red line) versus the optical detuning Ao, 
for unbalanced decay rates. Solid dots (•) are from the numerical integration of Eq. (Q with parameters Q,i — 5 x 10~'^r, 
7c = 5 lO^^F and ^2 = 5 x lO^^F. In (a) F32 = lO'^^F, F31 = F - F32. Note that the (SCJ33 shift is indistinguishable from the 
shift. In (b) F32 = IQ-^^F, F31 = F - F32. 



where 



S" = 1 



7.r 



3i 



5^,(^ = 1,2), 



2nj 

All =7ii (Ai +71A) , 
A22 = 722 (A2 - 72A) , 

^ _ »ir3i72Ai - nlr32iiA2 
7i72(f^ir3i + r!2r32) ■■ 



(22) 



and 



with 



7n = 7c7e// + 7ii (71 7e// - 7cAiA) , 

7I2 = 7c7e// + 722 (727e// + 7cA2A) , 



(23) 



7ii 



722 



21]hi + (A?+7?)r3i^ 

r31^2 + ^32^1 

21^172 + (Ai+7|)r32- 



(24) 



Fig. 2] reports the population resonance under various 
saturation conditions, matching the AT, EIT and FF 
hncshapcs of Fig. Notice that the EIT regime of 
Fig. mb) corresponds to the case of Rabi frequencies 
smaller than the natural decay rate of the excited state. 
The "Lamb-dip" like lineshape for the resulting quasi 
saturated transition can conveniently be observed in the 
three-level configuration. This method of spectroscopy 
without Doppler broadening was proposed and exper- 
imentally accomplished by Javan and Schlossberg in 
[gsI l69j . In such situation, the dip can be narrower than 
the homogeneous linewidth of the population resonance 
as in Fig. IDJb). 



The clock state populations depend strongly on the 
normalized branching ratio difi'erence T and on the Rabi 
frequencies driving atomic or molecular transitions. A 
numerical analysis of [T^I demonstrated the occurrence 
of a strong population transfer for unequal r3i, r32 decay 
rates. Fig.[S]shows the steady-state complete population 
transfer for fli ^ fl2 and T w 1 using the asymmetric 
decay rates associated to an alkaline-earth three-levels 
system as strontium atoms [ISjIlll ■ A large coherent pop- 
ulation transfer P22—P11 = ±1 is achieved when fli ^ 0,2 
for T - -M or 1^2 > Hi for T 1. 



B. Approximated frequency-shift of the 
two-photon resonance 

In an alkaline-earth frequency clock probing scheme, 
the large population transfer regime of Fig.[5]may be used 
to detect one lower state population, or the population 
difference between clock states. Thus, it is important 
to derive the two-photon shift also in this scheme. An 
analysis equivalent to the derivation of Eq. when 
Ai ^ A2 ~ Ao, we obtain, the frequency shift of pn. 
For P22 , we make a similar derivation also using the pop- 
ulation conservation condition. We obtain the following 
expression for (5ajii(Ao) and similarly for 61022(^0) 

(5a;ii(Ao) ^ Suj22i^o) ~ 

r? 



7?i 



All 



(A/-Aii)(7fiA/-Aiir2„) 



(r^//-7?i) 



2 ^2 



(25) 



For the specific radiative configuration of alkaline-earth 
species shown in Fig. [51 only the solution with the mi- 
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nus sign is needed but =F solutions generally refer to the 
extrema of a dispersive lineshape. Our standard choice 
of the laser detunings Ai = Aq and A2 = Aq — 5r 
introduces a very small difference in expressions of the 
frequency-shifts affecting each clock state population. 
For a quasi-resonant interaction when 7c7e/ / /^Iff ^ 1 
with Af/Teff 1 and Au/Teff <C 1, the expression 
may be simplified to yield: 



Su!ii{Ao) ^ Jw22(Ao) 



All 



7ii 

r 



2(1 



7f 



:) 



(26) 



The frequency shift versus the common mode optical de- 
tuning Aq affecting the pii(equivalently P22) resonance 
is plotted in Fig. HJa) and (b) for two particular ratios 
of the relaxation rates by spontaneous emission. In both 
cases, the shift of the two-photon resonance measured 
on the pii or P22 observables has a dispersive lineshape 
versus the optical detuning Aq. The slope is completely 
reversed owing to a nonlinear behavior when the ratio 
T 1—^ 1 as in Fig. [6Kb). A comparison with the frequency- 
shift of the excited state ^33 is also included in the figure. 
Notice the diflterence in the pi 1/^33 shifts for the case of 
a large asymmetry in the spontaneous decay rates. 



V. STEADY-STATE LINESHAPE OF RAMAN 
COHERENCE pi2 

A. The Raman coherence resonance 
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We are now focusing on the steady state Raman co- 
herence resonance pi2 = Re{pi2} + i /m{pi2} given by 



Re{pi2} = ri25^ 

/m{pi2} = Fis^^ 



(Ai - A2)A - 7e// 



(Ai-A2)[Ai-A2-A/] + r^^^ 



(Ai 



A-fc 



2 ' 
(27) 



(Ai-A2)[Ai-A2-A/] + r2^/ 



with 



ri2 



^^2^31 



32 



217iO 



lSi2 



(28) 



When the dipolc transition is allowed, the Raman coher- 
ence resonance can be detected in several manners. If we 
deal with alkaline atoms such as Cs or Rb, the hyperfine 
Raman coherence might be detected as a microwave emis- 
sion proportional to |pi2p, inserting the atomic medium 
into a micro- wave cavity [71[. Fig. [7] shows I/O12P versus 
the Raman detuning condition 6r when r3i = r32 = r/2 
at different values of Ag. The dispersive behavior of 
Rs{pi2} leads to a second resonant peak, appearing when 
Ao > 3r as seen in Fig. EKb) and Fig. EJc). 
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FIG. 7: (Color online) (a) Steady-state lineshapes of 
square modulus versus 5r using Eq. (|27p for Fsi = F 
fix = 5 X 10"^r, ^2 = 5 X 10~^r, and 7c = 0. In (a 
in (b) Ao = SF; in (c) Aq = 20r. Solid dots (•) are 
numerical integration of Eq. @. 



the|pi2p 

32 = r/2, 

Ao = 0; 
from the 
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FIG. 8: (Color online) (a) Frequency shift of the Raman coherence resonances observed on the |pi2p as derived from Eq. (|29|) 
(solid blue line) versus the optical detuning Aq. In (a) we have a symmetric radiative configuration with Fai = F32 = F/2 
where fii = 5 x lO^^F, ^2 = 5 x 10"^F and 7^ = 0. In (b) we have an asymmetric radiative configuration with F31 = F — F32, 
F32 = lO^^F where Hi = 5 x 10~^F, fl2 — 5 x 10~*F and ■jc — 5 10~''F. Solid dots (•) are from the numerical integration of 
Eq. dl). 



B. Approximated frequency-shifts of the Raman 
coherence resonance 



We derive here an accurate expression for the 
frequency-shifted resonance of |pi2p when Ai ^ A2 ~ 
Aq. a cubic equation is derived from the analytical dif- 
ferentiation of 1/512 1 ^ of Eq. (P7)) with respect to the Ra- 
man detuning Sr- Using Cardan's solutions, the three 
roots are written as: 



(5a;i2(Ao 
b 

3a 
with k = 



Cos 



1 



-Cos- 



3q 
P 



2k- 



(29) 



0, 1, 2, the p and q given by 
3ac - &2 



P ■■ 



3a2 
26^ - Qabc - 



27a^d 



(30) 



27a3 



and 



= -1- A 

= 3A (7e// 

= - 2 (7.^ A^ 
+ A/A(7, 



-leff), 



(31) 



'^=re//A(7c + 7e//) + A/ (7c A 



Eq. allows us to obtain the frequency shift as a 
function of the common mode optical detuning Aq plot- 
ted in Fig. [HI However an estimate of that shift is ob- 
tained looking only at the real part of the coherence solu- 
tion which mainly describes the lineshape emission. From 



the square modulus of the real part |_Re{pi2}p, simple 
cubic solutions for the Raman coherence frequency-shift 
can be derived as: 



or 



(32) 




SuJi2{Aq) 



leff 



(33) 



When ATeff/jeff < 1 and A//Fe// <C 1, the Ra- 
man coherence frequency-shift expression corresponding 
to the maximum emission becomes: 



Af 



(34) 



For that case, we recover the usual dispersive shape re- 
lated to the light-shift Als affecting clock states. The 
frequency shift versus the common mode optical detuning 
Aq derived from Eq. ((29)) is shown in Fig. |8l for the case 
of a symmetric radiative configuration with F31 = F32 
while in (b) for F31 ^ F32. The central dispersive fea- 
ture, related to the Raman shift expression Af, corre- 
sponds to the maximum of the coherent emission. Other 
branches of the shift correspond to the extrcma of the 
second resonance appearing for Aq > 3F as from the 
lineshape simulation of Fig. El^b). A direct comparison 
of Fig. [HUb) with frequency shifts reported in Fig. [6] with 
similar conditions, yield to Raman coherence shifts larger 
by more than an order of magnitude than population fre- 
quency shifts. 
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VI. DARK RESONANCE FRINGES 



A. Pulsed regime lineshape 

The clock operation may be based on a pulsed Raman- 
Ramsey scheme, illustrated in Fig. |9l with beating oscil- 
lations observed whichever variable is monitored. This 
detection approach originally introduced in [53l. was re- 
fined in [Ho, Im and discussed in refs (5345a | . It allows 
to reach a higher precision in the clock frequency mea- 
surement, as typical of the Ramsey fringes. The present 
work focuses on the laser pulse scheme where the first 
pulse is long enough to allow the atomic or molecular 
preparation into the dark state superposition and the 
second delayed short pulse probing that superposition. 
Because the Tp pumping time of Eq. is required to 
reach the steady-state atomic or molecular preparation 
into the dark state, the length of the first pulse should 
satisfy r ^ Tp(Ao) ;» r3j^^,r32^. From a mathemat- 
ical point of view the t i— > oo steady state solution 
of the three-level optical Bloch equations may thus be 
used as initial condition for determining the evolution 
at later times. At time r the laser fields are switched 
off in order to allow for a free evolution over the time 
T. Finally a readout operation is performed by apply- 
ing a short pulse whose time duration r™ is limited by 
^31^^32 ^ Tpi^o)- In this adiabatic regime 

where dp33{t)/dt = dpi3{t)/dt = dp2-i{t)/dt = 0, the 
short probe pulse operation is well described using the 
T,„ t-^. limit. For a readout pulse as long as the prepa- 
ration pulse, all atoms or molecules are repumped into a 
new dark state erasing interference fringes. Interferences 
fringes are detectable on all observables as a function of 
the T time delay, with very short readout pulses required 
for optical coherences and the excited state population 
fraction. Instead longer probing times are required for 
monitoring fringes on lower state populations due to slow 
time evolution of the clock states and the Raman coher- 
ence. 

An analytical expression of the popula- 
tion fraction in the excited state can be established 
from Eq. ^ in the asymptotic limit of a long 
preparation pulse and a vanishing read out pulse, 
i.e. at the end of the free evolution time. In 
this regime, a pulsed DR expression p^^iT) and re- 
lated transmission parts of the /m{/9i3}(T), Im{p23}{T) 
optical coherences can be rewritten in the exact 
ay [l-H/3„ |pi2|e-'>''='^cos((5rT-$,j)] form: 



P33(T) = a33X (1 +/333|pi2|e '^'^'^COs(J^T-$33)) 



"33 



J33 



»i(An + l)/r72 - »f (An - l)/r7i 

i-ff(r!?/7i + f}2/^2) 



(35) 



^733 + MA33 



X » T 

Preparation time 



Ramsey time Probmg time 

T„ « T„ T„ » t:„ 



Dark Resonance Dark Resonance fringes Dark Resonance 

,i. 





FIG. 9: (Color online) Pulsed Dark Resonance detection 
scheme to perform high resolution spectroscopy of three-level 
A systems. T is the Ramsey time when both laser light 
fields 0,1 , 0,2 are switched off. The first pulse is long enough 
T ^ rp(Ao) to reach the stationary regime. During the sec- 
ond pulse, the probing time can be short r,,, <C rp(Ao) to 
observe Dark Resonance fringes or long ^ Tp (Ao) to re- 
cover a cw Dark Resonance as a new preparation stage for 
the next pulse. 



/m{pi3}(T) =ai3 X (1 + /?i3|pi2|e-'^=^ cos ((5,T - $13)) 



ai3 =^ (3q;33 + (An - 1)) , 
271 



(36) 



^713 + A^A13 



J13 



ai3 



/m{p23}(T) = a23 X (l + /323|pi2|e-^=^cos(<5.r-$23)) , 



V12 

"23 =77^ (3a33 - (An -I- 1)) , 



272 



(37) 



/3: 



M723 + A'A23 



23 



^23 



where An = P22 — pii is the steady-state clock state pop- 
ulation difference obtained from Eq. (j2ip and \pi2\ given 
by Eq. (|27p . The following quantities were introduced in 
the above expressions: 



M733 



MA33 



2VtiVl2 f (71 +72)/7i72 

r 1^1+ 3(l]2/^^+f^2/~) 
2VI1VL2 ( Ai/7171 - 1^2/1212 

l + f(0?/7i+f^i/72) 



(38) 



2,9,1 

2 tT 
317i 



^2 

2 71 71 



033 



o ^1 

/^A13 = 77 — /^A33 - —: 

2 71 7171 



(39) 



12 



2.0x10"" 



O 1.5x10 



•2 LOxlO"" 



5.0x10' 




-2.0x10 



-1.0x10* 00 1.0x10" 

Raman detuning 8 IT 



2.0x10 



0.15 




-4-3-2-10 1 2 3 4 
Raman detuning 5 IT 



2.0x10 



o 1.5x10* 



■2 1.0x10 



5.0x10 




-2.0x10"" -1.0x10"" 00 1.0x10"" 
Raman detuning 5/r 



2.0x10 



FIG. 10: (Color online) (a) Dark Resonance fringes in the 
weak field limit, (b) Frequency span of the lineshape. Lines 
from Eq. (|35[) and dots (•) from Bloch's Eqs. Common 
parameters: Fgi = F32 — F/2, Ao = 0, 7c = and T — 
lOrp(O). In (a) and (b) Rabi frequencies fli = fl2 = O.OOSF, 
free evolution time T = 10 Tp(0), probe time ~ 0.004 Tp(0). 
Very good agreement between Eq. Q and Eq. (|35p results. 



3^2 



Hi 



M723 — 77-/^733 — 
2 72 ' 72 

MA23 — 77— A*A33 + "1 —■ 

2 72 7272 



(40) 



In the non asymptotic limit where r T„i ^ Tp, 
stationary solutions of \pi2\ and An have to be re- 
placed by their transient expressions \pi2{T,Tm)\ and 
An(T, Tm). When, T M> 0, the lineshape expression 
from Eq. ([55)1 is formally equivalent to Eq. ([5]). A 
pulsed DR lineshape originated is plotted in Fig. [TOla) 
for small values of Rabi frequencies with a magnified 
span on the central fringe in Fig. [TUT b). In these plots, 
the agreement between Eq. ([55]) (solid red line) and 
the Bloch's equations (dots •) is very accurate. We 
have also studied the role of the Ramsey time T on 
the DR lineshape and amplitude. Fig. [TlT a)-(c) shows 
the signal amplitude versus the Raman detuning for 
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FIG. 11: (Color online) DR lineshapes in strong laser fields 
computed from Eq. (|35|) versus Raman detuning for different 
Ramsey time T. Parameters i}i = — 0.2F, T31 = F32 = 
F/2, Ao = and 7^ = 0. Solid dots (•) from Bloch's Eqs. © 
with r = 10 rp(0). In (b) and (c) r„ = 0.75 rp(0). 



different values of T. The fringe amplitude is always 
twice the amplitude in the cw regime when 2T > ^^ff 
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FIG. 12: (Color online) Asymmetric fringes observed on the imaginary parts of optical coherences, /7Ti{pi3}(r) in (a) and 
Im{p2s}{T) in (b), and on the excited population pss, in (c), versus Raman frequency detuning. In (d), (e), and (f) expanded 
views of the central resonance fringe for each observable. Parameters Fsi = O.SF, F32 = 0.2F, Ao = F, 7c = 0, Q,i = 5.10~^F, 
512 = 2.5 X 10~^F, free evolution time T = 10 rp(0), r > lOrp(F), ~ 5.5 x 10~^rp(F). Solid lines from the analytical forms 
of Eqs. (|36|) . (|37p . and (|35p . respectively, in good agreement with the solid dots (•) results from Bloch's Eqs. Q. 



and 7c = 0. A careful comparison between Fig. 1101 and 
Fig. fTTlb) and Fig. [TTTc) show that for Rabi frequencies 
comparable to the F3i,r32 relaxation rates, i.e., when 
saturation becomes important, the limit Tm '-^ used to 
obtain Eq. ((35)) is no longer valid and the effect of the 
readout duration pulse length has to be included into 
the analysis leading to a slight reduction of the fringe 
amplitude. 



B. Stationary atomic/molecular Raman phase 

The fringe signals appearing in the atomic observables 
arc produced by the cos{6rT — terms appearing in 
the above expressions. Within a narrow region around 
the resonance, the lineshape is entirely described by a 



phase given by 



33 



13 



Arctan 



Arctan 



^^733 Im{pi2} - AtA33 Re{pi2} 



M733 Re{pi2} - 
^^13 Im{pi2} 



/XA33 Im{pi2} 

- PA13 Re{pi2} 



$23 = Arctan 



M713 i?e{pi2} + /XA13 Im{pi2} 
p^23 Im{pi2} - PA23 Rejpu} 

M723 Re{pi2} + PA23 Im{pi2} 



(41) 



In the adiabatic regime, for a first pulse producing 
a steady state and a second pulse with a short dura- 
tion, the IP12I steady state solution of Eq. (|27)) leads to 
cw expressions of $y (i,j = 1,2,3) where all time de- 
pendances are absent. For clock engineering, it is im- 
portant to quantify the shift of the central fringe when 
the common mode optical detuning Ai ^ A2 Ag is 
scanned around the Sr = resonant value. Within the 
limit of a long first pulse and a short second pulse, in 
the weak field regime and for a sufficiently large Ram- 
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FIG. 13: (Color online) Exact numerical tracking of fringe 
frequency-shifts Suii^, SUJ23 and ^0^33 versus Ao derived from 
Eqs. (I35j, dMJ and dSZll. Parameters Tsi = O.SF, r32 = 0.2r, 

7c = 5 X lo^-^r, fii = 5 X lo^^r, n2 = 3.7515 x lo^^r, t = 

10 Tp(0). Solid dots (•) from Bloch's equations Q where we 
use T > lOrp(r), Tm ~ 5.5 x 10~'^Tp(r). Discrepancy between 
solid lines and dots are here due to the non vanishing Tm 
readout time required with a numerical integration of Bloch's 
equations. 



FIG. 14: (Color online) Exact numerical tracking of the Dark 
Resonance fringe frequency-shift 5uj% versus Ao derived from 
Eqs. ((351). Solid dots (•) are plotted from Bloch's Eqs. (g]) 
with parameters as in Fig. 1131 



sey time (T ^ rp(Ao),T„i), the central fringe frequency- 
shift Sojfj should be connected to the $jj phase accumula- 
tion. Time-dependent Raman frequency-shifts produced 
by short preparation pulses were examined in details in 
ref. [tI, [t^ using the Bloch vector model, but only for 
the excited state population. A complex wave-function 
formalism was also proposed for short pulses in (H^j and 
extended in ref. [s^ to derive analytical time-dependent 
expressions of any frequency-shift affecting the central 
fringe related to each possible observable. 

Fig. [IHa)-(c) reports the fringes appearing on the 
imaginary parts /m{pi3}(T), Im{p23}(T) and on the ex- 
cited state population p^^iT) for an asymmetrical radia- 
tive configuration with T^i = Q.8T and = 0.2r, a de- 
tuned laser excitation (Aq = F) with no Raman decoher- 
ence {jc = 0). Plots from Fig. [T^d)-(f) report a magni- 
fied span on the central fringe. The plots in Fig.[T^a).(b) 
for the optical coherences show asymmetrical lineshapes 
with the central fringe being blue or red shifted from the 
exact Raman resonance, a behaviour different from that 
observed in the cw regime. For the excited state fraction 
plotted in Fig. \W[c) under same conditions, the central 
fringe is not frequency shifted, as in the cw regime. These 
results point out, for the first time, that the imaginary 
parts of optical coherences, when individually probed, 
have a hneshape different from that observed on the ex- 
cited state response due to asymmetric decay rates by 
spontaneous emission. 

This different behaviour is also confirmed by the plots 
of Fig. [13] focusing on Dark Resonance fringe frequency- 



shifts versus the common mode optical detuning Aq when 
assuming a small decohcrence term 7c = 5 x 10~*r. 
Fig. US] are numerical tracking of frequency-shifts from 
the Eq. (j35p analytical form compared to the integration 
of Bloch's equations (solid dots) for a particular ratio 
between Rabi frequencies. Oscillations around the nu- 
merical track of frequency shifts are observed for large 
optical detunings due to the non vanishing Tm readout 
time used with Eqs. (|4|). Fig. [T4l evidences a very weak 
slope near optical resonance; such dependence does not 
appear on the imaginary parts of the optical transmis- 
sion. In absence of the jc decoherence the excited state 
dependence was recently discussed in works [t^, [ill deal- 
ing with the time-dependa nt p art of the Raman- Ramsey 
fringes reported earlier in [72l [73|. 

A perturbative expansion of Eq. p5|) in the Raman 
detuning parameter is required to derive the correct non 
linear behavior of the frequency-shift versus the Aq de- 
tuning, as reported in Eq. (jB9p of Appendix B. Those de- 
pendencies (solid dots) are plotted in Fig.[T5]for ^12/^1 ^ 
0.75 and compared to the numerical tracking shift of the 
fringe minimum from Eq. (j35p (solid lines). A non lin- 
ear behaviour with small oscillations around the optical 
resonance is correctly described. Increasing the T time, 
there is a small rotation of the frequency-shift around the 
Ao = common mode detuning as shown in Fig. [TSlb). 
This dependence appears only for a non vanishing 7^ de- 
coherence rate. For a given value of T, at specific val- 
ues of the optical detuning, an exact cancelation of the 
frequency-shift takes place as seen in Fig. [TSlb). 
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VII. CONCLUSION 

Three-level systems interacting with two coherent 
laser fields give rise to many phenomena, such as 
Autlcr-Townes doublet, Dark/EIT Resonance, and Fano- 
Feslibach transition. In all of them the Raman coherence, 
playing an essential role, is highly sensitive to the param- 
eters of the laser interaction and of the decoherencc pro- 
cesses. The present work explored how those parameters 
may be tuned for future optical clock devices providing 
sensitive measurements of energy levels shifts in neutral 
atoms clocks based on either fermionic or bosonic atomic 
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FIG. 15: (Color online) (a) Expanded view of the fringe 
frequency-shift numerical tracking from Fig. [T?] for var- 
ious values of the small decoherence 7c term, (b) Exact Nu- 
merical frequency-shift 5a;|3(Ao) with different Ramsey time 
T for the particular ratio S72/f^i ~ 0.75 where f2i = 5 x 10~^r 
and 7c = 5 X 10~*r. Solid dots in (a) and (b) cases are plot- 
ted from the analytical first order expansion of the Raman 
shift using Eq. HB9|I . Note the oscillating pattern with small 
amplitude. 



species for example in dipolar traps |76| or in trapped ion 
clocks [77j . 

The three-level phenomena are well described by the 
formalism of the Bloch's equations in a semi-classical 
density matrix representation. From the exact resolu- 
tion of Bloch's equations we have derived general an- 
alytic expressions of the resonance lineshapes observed 
in the steady-state of different atomic observables. We 
have examined the linewidth and power broadening of 
the two-photon resonance. The precise dependence of 
the frequency shift associated to the Fano-Feshbach tran- 
sition or to the FIT resonance, not appearing in a per- 
turbation treatment, was discussed. We have exam- 
ined the resonance fluorescence and frequency-shift for a 
pulsed laser configuration in the adiabatic regime. The 
analytical but asymptotic solution allows us to write 
the lineshape solution of the narrow quantum resonance 
leading to formation of Dark Resonance fringes. The 
pulsed sequence overcomes the power broadening mech- 
anism of the continuous-wave excitation while allowing 
high contrasted signals in a saturation regime. In the 
weak field limit, signals in pulsed regimes are two times 
the continuous-wave signals, except for the decrease due 
to the Raman coherence relaxation within the atomic free 
evolution. The lower limit to the resonance linewidth is 
1/2T to be compared to the Fe// limit in the steady-state 
case. The resonance shift is still proportional to the 7c 
decoherence as for the steady state regime, but it is now 
diluted over the Ramsey time. 

An important and original result of our analysis is that 
different atomic or molecular quantum observables (ex- 
cited population, clock-state populations and Raman co- 
herence) experience different non-linear lineshapes and 
therefore different shifts of the clock resonance. Depend- 
ing on the atomic parameters, the cw shifts of those ob- 
servables may greatly differ in slopes, magnitude or line- 
shape. The shift amplitudes are also strongly related to 
the values of the decay channels, either balanced or un- 
balanced. For instance, while a large clock-population 
inversion is produced for unbalanced decay channels, a 
dispersive lineshape is associated to its frequency shift. 
The detection of the Raman coherence exhibits a larger 
contrast, but it suffers usually of a systematic shift larger 
by an order of magnitude than that associated to the 
population detection. In the case of a very small decay 
rate for the Raman coherence, destructive interferences 
strongly reduce the cw frequency-shift sensitivity of the 
excited state compared to others observables. Thus the 
choice of the observable is very important for the proper 
operation of a three-level atomic or molecular clock. Ad- 
ditional technical constraints are associated to the re- 
quested detection tools of the chosen observable. For in- 
stance tracking the excited state population fraction by 
monitoring spontaneous emission or light transmission 
represents a sort of non-demolition quantum measure- 
ment avoiding the destructive readout associated to the 
lower-state clock projection. 

The present solution can be extended to the case of 
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FIG. 16: (Color online) Definition of the Wij transition rates 
within the three-level A configuration 



a train of laser pulses having different frequencies, op- 
tical detunings and phase steps as in refs [zl, [T^I, in 
order to design more elaborated combinations of optical 
transient nutations and free evolutions and to explore 
more efficient detection scheme for very narrow transi- 
tions. For a treatment of the atom or molecule motion 
in matter wave interfcrometry [sof , for instance based on 
a pulsed EIT/Raman interaction with a resonant excita- 
tion scheme and using stimulated photon recoils to create 
a beam splitter, the Bloch's set of Eq. ^ requires both 
the introduction of the recoil shift into the laser detunings 
and the recoil spread among different atomic momentum 
classes produced by spontaneous emission yielding to a 
full set of quantum coupled equations [sij . 
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From the steady-state solutions of rcfs. [i^li^l we obtain 

W12W23 + Wi3{wi2 + W23) 



P33 = 
P22 = 
Pll = 



3wi3Wi2 + r32W"ll -I- r3iW22 + 3W23W33 ' 
W12W13 + r3iWi2 -I- r32(wi3 -I- W12) + ^23^33 

3wi3Wi2 + r32Wii -I- r3iW22 + 31023^33 
W12W13 + r32l«i2 + r3i(lt;23 + W12) + W23W33 



3wi3Wi2 + r32Wil + r3ilV22 + 3^23^33 



(A2) 



where the Wij transition rates between each pair of states 
are defined in Fig. 16 and where 



W22 
W33 



Wi3 + 2wi2, 
W23 + 2wi2, 
Wi3 + Wi2. 



(A3) 



The expressions for the two-photon {wi2 cx fli^2) ^-^^^ 
one-photon (1013 oc fl1^w23 oc ^2) transitions rates are 



W12 =2Re{ 



^31'?2lC3 



Wi3 =2Re{- 



^l£,3l + ^2^32 
— } - -^12, 
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W23 =2i?e{- 



521 + 



(A4) 
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Appendix A: Generalized multi-photon steady-state 
rate solutions 



The detuning parameters are complex numbers de- 
pending on the system parameters 



61 = +7i, 

62 = -iA2 + 72, 

^21 = «A2i + 7c = «(Ai - A2) -I- 7c 



(A5) 



The above expressions of the one and two-photon rates 
point out the light-shift contributions to the eigenfre- 
quencies of the three-level A system. 



The lineshape expressions of Eq. ([8]) and Eq. ((2T|) may 
be recast in a different mathematical form using the 
generalized multi-photon rate solution of ref. [1^, later 
employed in ref. [43| to establish conditions for the co- 
herent population trapping in the steady-state regime. 
That analysis is based on three coupled equations for 
the Pii{Z), with i = 1,2,3, occupation probabilities of 
the three levels, equations written in the Laplace space 



Appendix B: Dark Resonance fringe frequency-shift 
1. First order expression of Suj^^ 

Eq. psp can be recast in the following form: 

P33(r) = a33 + S33COs(C/), (Bl) 
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where 



B33 = M|Pi2|e '''='^, 

/U = "33/333, 

C/ = (5,r - $33- 



(B2) 



To establish the frequency-shift, a track of the extremum 
of Eq. ([55]) produced by a differentiation versus Sr leads 
to the following expression: 



P33 



^33 



^Bg + Bl^U'Hin{U + 633), 



with 



(B3) 



If we neglect 633(0) and 7733(0) terms, the discrepancy is 
respectively +6% and —13%. If we neglect the 633(0) or 
the 7733(0) term, the error is below +0.4%. Thus at the 
1% level accuracy and at a common mode detuning small 
compared to spontaneous emission rates, we can rewrite 
the frequency- shift expression as 



^33(0)- 633(0) -7?33(0) 

r-$33(0) 



(B9) 



733 



733 



= - arctan(^f-), if T - '^'^^ < 0, 



= — arctan(' 



B. 



(B4) 
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'B33U' 

The P33 derivative vanishes at 

U + 633 + 7]33 = 

where 



TT, if T- $33 > 0. 



?733 = arcsin( 



a 
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Rewriting the cancelation condition as: 

SrT - $33 + 633 + ?733 = 0, 



(B5) 



(B6) 



(B7) 



Derivative terms 



This subsection reports the functions appearing in 
Eq. (jB9[) and required to determine the frequency-shift 
Jwgg as a function of the common mode detuning Ag 
and the system parameters. The additional functions 
may be obtained from the first order derivatives of the 
/i, <I>33, IP12I, a33 functions calculated at Sr — 0. In or- 
der to simplify the mathematics, those functions are 
here reported for the case of a pure radiative process 
7i = 72,71 = 72, and supposing Ai = A2 = Aq. 



^^•33(0) = —Arctan 



7cA 



lleff 



(BIO) 



we apply a first order expansion of Eq. (jB7[) into the Sr 
Raman detuning leading to the following Suj^^ frequency 
shift: 



(JWgg — Sr ~ 



^33(0) -633(0) -7733(0) 

r -<i>:,3(0) + 6^3(0) + 77^3(0) • 



(B8) 



Aq (fir + snl - 



(Bll) 



^33 



(0) = 



AApnl (rfi + An (rfi + enj)) 



71 



Aii^" + A2252' - 7cA722'5'' + IcAo—^ (711^" + 722-5'') 



(7Vr + 37i(i7?-ff^i))r,2ff 
2nl 



4AoSnl 



1 + 2S 



r2 



2 f r, , 7ir32 



4Aoj2^S^^T32 (2AoA7e + T-f.s) 

r32 (4A2 + r2) + 4rnl 



2 I r,, 7~ir3i 



An A 



8Aon?(r2^-7,^)' 



y ^ 4A2 + (1 + 2S')r2 

2Ao7|2r32^'' 



ml 



(B12) 



*33(0) 



A%c (27?7i - 2Ag7c + 7171 7c) + leff (27171(71 + nljc) - 2A^7e// + 7i7i7e// 



27?7i (■ 



A\2 



(B13) 
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|P12| (0) = |pi2|(0) 



8Ao (4A^7,5 + r (7^^^ + 37e7eff>g - + 25))) 

(i + 25 + 4A2/r2)'r3r2jj 



7off + 2r!i72/r 



(A272 



I — 



BAof^i 



r2(i + 2S' + 4A2/r2) 



(B14) 



These perturbative expressions were used to determine 
the (5u;33 frequency-shift plotted in Fig. [TSl when Aq < F. 
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